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Abstract 

The problem of defining a gauge invariant effective potential with a strict 
energetic interpretation is examined in the context of spontaneously broken 
gauge theories. It is shown that such a potential can be defined in terms of a 
composite gauge invariant order parameter in physical gauges. This effective 
potential is computed through one loop order in a model with scalars and 
fermions coupled to an abelian gauge theory, which serves as a simple model of 
the situation in electroweak theory, where vacuum stability arguments based 
on the scalar effective potential have been used to place lower bounds on the 
Higgs mass. 
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1 Introduction 



Effective potential calculations in electroweak theory have been the subject of sus- 
tained interest as a consequence of the observation that phenomenologically interest- 
ing lower bounds on the Standard Model Higgs mass can be obtained from vacuum 
stability considerations (see [|l| for a comprehensive review of work through 1989, 
for a review of more recent work). Although the gauge interactions do not play 
a qualitatively significant role in these estimates, it is important to include them if 
one wishes to obtain as precise a bound as possible. Unfortunately, the conventional 
effective potential is an inescapably gauge-variant object, as the scalar field must 
transform nontrivially under the gauge group in order to break the symmetry spon- 
taneously in the first place. While general theorems ensure that the value of the 
effective potential at local extrema is a gauge-invariant quantity, the location of the 
extrema and the behavior of the potential between extrema can vary widely from one 
gauge to another. This makes it difficult to reach unambiguous conclusions on the 
basis of vacuum stability arguments 0]. 

There is a long history of attempts to formulate spontaneous symmetry breaking 
(SSB) in terms of gauge-invariant order parameters 0, 0]. In the case of dynamical 
symmetry breaking, where the spontaneous breakdown arises from the appearance 
of a vacuum expectation value for a gauge-invariant operator, such order parameters 



have necessarily been taken to be composite operators p, [T^. Thus one is led to 
the introduction of an effective action for composite gauge-invariant fields, defined 
in the usual way as the Legendre transform of the connected generating functional 
for n-point insertions of the composite operator . Such approaches can further be 
divided into those based on a local composite operator, and those where the order 
parameter involves a bilocal operator, i.e. a product of fields at separated space-time 
points. In the former case it is well-known [^, [1^ that the additional subtractions 
needed to render the n-point functions of the local composite operator finite vitiate 
the energetic interpretation of the effective potential (obtained as the translationally 
invariant limit of the effective action). This is obviously fatal for any attempt to 
study vacuum stability using such potentials. Other authors have employed bilocal 
operators |^ with gauge-invariantizing factors, but in non-physical covariant gauges 
where a rigorous energy interpretation is again lacking. 

In this paper we present a complete calculation to the one-loop level of a com- 
posite effective potential based on a gauge-invariant order parameter in a physical 
gauge (Coulomb) where the energetic interpretation of the potential is preserved. 
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The bilocal operators used involve a smearing function which is largely arbitrary. 
The smearing dependence of the results can however be understood in a completely 
physical way, and is found to be numerically insignificant in the regime of interest. 
The formalism needed is illustrated first in a simple scalar model (Section 2), and then 
extended to the Higgs-abelian model coupled to fermions (Section 3). The effective 
potential computed in Section 3 is renormalized and the final explicit finite result 
for arbitrary smearing functions given in Section 4. Section 5 contains an explicit 
analytic evaluation for the special case of a smearing function with a sharp cutoff in 
momentum space. In Section 6 we apply the composite operator formalism to the 
issue of electroweak vacuum stability bounds. Renormalization group improvement 
of the composite effective potential is also discussed here. Some explicit numerical 
results are presented in Section 7. The issue of extension of the formalism to non- 
abelian gauge groups leads naturally to the subject of composite effective potentials 
in axial gauge. The choice of an appropriate bilocal operator in this case involves 
a number of subtleties at both the perturbative and nonperturbative level. Some of 
these issues are discussed in Section 8, although a full computation at the one loop 
level is deferred to a future publication [|T3| . 

2 One-loop effective potential for a composite op- 
erator: a simple example 

In order to establish notation, and to remind the reader of some features of the loop 
expansion of effective potentials which arise when the order parameter is a bilocal 
composite field, we first consider a simple model of a single real scalar field with 
Lagrangian 

£ = U.ct^d^^ct^ + ^m^02 _ p(0) (1) 

Conventionally one studies symmetry breaking in this model by constructing an ef- 
fective potential defined |T^ as the Legendre transform of W{j), the generating func- 



tional of connected graphs: 

W{]) = -i log Z(j) = -i log y V0 /(^-^■^)'^'^ (2) 

The corresponding Legendre transform, r(0) (where now represents a classical field) 
has a straightforward interpretation in terms of one particle irreducible (IPI) graphs 
and is easily evaluated graphically at one-loop order. If we wish instead to add a 
source for a composite bilocal operator, the graphical interpretation is more compli- 
cated, and it turns out to be easier to construct the one-loop potential by a direct 
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semi-classical expansion. The calculation is also readily performed for a composite 
operator in which the field point splitting is smeared spatially (in all directions) with 
a function K{f), which is essentially arbitrary, except for (a) being nonsingular at 
the origin r = (in particular, there must be no delta function singularities there) , 
and (b) having a positive Fourier transform K{p). The reason for the first condi- 
tion is well-known- ri-point functions of local composite fields will require additional 
subtractions which correspond to nonlinear source terms which ruin [|I^ the energy 
interpretation of the effective potential. The positivity requirement on K{p) ensures 
that if we define our composite field as 

X{x) = J (j){x + r/2)K{r)(p{x - f/2)rf¥ (3) 

then expectation values of x ^"^^ finite (once the field is renormalized) and positive. 
The latter property is easily seen by inserting a complete set of states between the two 
fields in < 0|x(a;) |0 > and using translational invariance. This is in distinction to the 
case of a local composite field, where the subtractions needed to define the composite 
operator would ruin the positivity, as well as introducing the need for counterterms 
nonlinear in the source which ruin the energy interpretation of the effective potential. 
Here we are assured that the appropriate domain for the effective potential V{x) is 
X > 0. Finally, it will be convenient to normalize the smearing function K{r) by 
/ dfK{f) = 1. This ensures that at the classical level, x reduces simply to </)^ in the 
translationally invariant case 4>{x) = const. In momentum space, this means that we 
take Ki^p = 0) = 1. A convenient choice is the Gaussian K{p) = e~T^I^'^P^^ where 
the parameter p is roughly the inverse smearing separation of the fields in coordinate 
space. However, it is not essential that i^' be a smooth function- one may also use a 
sharp cutoff in momentum space, with K[p) = 6{p — \p\), which has the advantage 
that the one-loop integrals for the composite effective potential can be performed 
analytically. As we shall see below, the qualitative results are similar in both cases 
when p takes a physically sensible value. 

The field x is not an order parameter for symmetry breaking in the conventional 
sense, vanishing exactly in the symmetric phase and giving a non-zero expectation 
value in the broken phase. Evidently, x 7^ even in the symmetric phase. However, 
we shall see below that if the smearing scale ^ is taken large compared to the inverse 
mass gap in the theory, the value of x the minimum of the potential does go 
to zero in the symmetric phase, while approaching the square of the conventional 
vacuum expectation value of in the broken phase. Consequently, the scale (if any) 
at which symmetry breaking occurs can just as well be studied by looking for global 
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minima of the potential V{x), provided that such a potential is defined in a manner 
consistent with an energy interpretation. This means that composite operators must 
be point-split (to avoid counterterms nonlinear in the source) and the point-splitting 
must be spatial only, thereby maintaining the Schrodinger picture derivation of the 
energetics of the system perturbed by a source jx- 

The calculation of the Legendre transform T{x) of W{j) 

W{j) = -i log Z{j) = -i log J V(j)ei Iic-3x)d^^ (4) 

through one loop is most easily performed by saddle-point techniques. As we are in- 
terested in the result only through 0{h), for the calculation of W{j) we need only keep 
Gaussian fluctuations around the saddle-point 0o of the source-augmented action: 

(□ + m2)0o(x) + P'(0o) + 2 1 dfj{x - ^)K{r}Mx - f ) = (5) 

Passing to the Legendre transform will require the elimination of the source j{x) in 
favor of the transform variable x- Taking the Gaussian smearing 
be specific, one may write formally (by Taylor expanding the fields 0(x ± |r*) in the 
definition (^) 

x(x)=0(x)2--L|V</)(x)|2 + ... (6) 

Evidently, as asserted previously, in the translationally invariant limit we recover 
simply X = 0^- Taking P(0) = Ao0^/4, the saddle-point equation then implies 

jo = -^{ml + \o<i)l) = -^{ml + Xoxo) (7) 

in the translationally invariant limit. Here, the subscript "0" on the fields 0, x refers 
to the leading order in the h loop expansion, while the subscripts on m, A remind us 
that these parameters are bare ones, with mo = mn + h ■ (1 — loop counterterms). 
The leading contribution to the effective potential (defined as the effective action, 
or Legendre transform of W{j) per unit volume in the translationally invariant limit 
where x{^) = constant), is obtained by evaluating the classical action at the saddle 
point 00 and eliminating 0o in favor of x = 0o + 0{h). As the saddle-point is by 
definition an extremum, the 0{h) correction in x only affects the result at 0{h'^) and 
so may be neglected to 1-loop order. Not surprisingly, we recover at order the 
expected tree potential, this time written as a function of the classical composite field 
X 

V,Ux) = lrnlx+^x' (8) 
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= ^ / 7^1(1^ + ml{l - Km + 3AoX - AoX^(p)) - In (p^)) (11) 



The one-loop contribution is obtained by integrating out the quadratic fluctuations 
around the saddle-point field 0o- Unlike the usual case where the source is coupled to 
the elementary scalar, so that the source term does not contribute to the quadratic 
part, here the source term acts as a momentum-dependent mass term at the one-loop 
level. Defining Ptot(0) = + JX the one-loop integral gives a contribution to the 
effective potential 

V^i-ioop = -^Tr In {U + ml + (9) 

For translationally invariant saddle points, the determinant in (P) reduces to (we are 
still in Minkowski space) 

Tr In (□ + + i^^J = j ^Mp^'^I' ^"(<^o) - 2joi^(p)) (10) 

where jo is to be written in terms of x using (0), and we may replace Xo by x 
everywhere in (^) as we are already at 0{h). After a Wick rotation to Euclidean 
space, the 1-loop contribution may now be written 

1 f d^p 
(2^ 

The second term subtracts off a field-independent divergence which is physically in- 
significant. Although we continue to write the 1-loop contribution in terms of the 
bare parameters tuq, Aq, these may be replaced with renormalized parameters in the 
1-loop term, which is already of order h. It will be convenient to perform the renor- 
malization of the theory using the MS scheme, as we shall shortly be moving on to 
gauge theories. In this scheme the relation between bare and renormalized quantities 
is (to order h) in 4 — e dimensions 

Ao = A^ + ^^A^A, (12) 

ml = mlil + h^XnA,) (13) 
2 

A, = 7 + ln47r (14) 

e 

where 7 is the Euler constant and the factor of h makes the 1-loop order at which we 
are working explicit. 

As the smearing function only involves the spatial components of the four vector 
p, we first integrate out the energy component explicitly, leaving a 3-e dimensionally 
regulated integrand: 

^1-loop = lhj ^^S^/f + 2\nx + {ml + A«x)(l - K{p)) - \p\) (15) 
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where we can replace bare by renormalized quantities as we are already at 1-loop 
order. The divergent parts of this expression may be isolated by writing 



Vl-loop = V^ri + Kub (16) 

^fin = \n J ^(Vp" + 2Ai?x + {ml + XrxKi - km - \p\ 

^ (m^ + 3AfiX) + ^7^^A^(^fl + 3Ai?x)') (17) 



- le^'i (27r)3-^ (p2 + ^2)3/2 

Here is an infrared cutoff introduced to ensure that T^in remains both infrared 
and ultraviolet finite when e is taken to zero. The first subtraction term vanishes 
in dimensional regularization. The subtraction term is easily computed in 3 — e 
dimensions: ^ 

l^sub = -T^Kml + 3A^x)'(A. - In \) (19) 

with /i the usual renormalization scale which is introduced by the dimensional contin- 
uation. Adding together the tree contribution (|]) (with the bare quantities replaced 
by the renormalized ones using (|T2|)) one finds the renormalized composite field po- 
tential through 1-loop, up to an irrelevant additive constant, to be 



Kcn(x) = \mlx + -MX^ + VUX) + ^("^R + 3ArX)' In ^ (20) 



The logarithmic dependence on rj in the last term is cancelled by a similar dependence 
in l^n, so we can take rj to zero after the integrals are performed. 

The local limit for the composite field corresponds to taking p in K{p) = e"^/^^''^) 
to infinity, i.e. K ^ 1. If we do this, a quadratic divergence reappears in Vfin, as 
the x-dependence at large momentum no longer matches that of the counterterms 
needed to renormalize the theory. This is a well-known difficulty, much discussed 
in the literature. Alternately, if we take p small, corresponding to separating the 
field points in the bilocal composite by large distances, the potential simply goes over 
smoothly to the effective potential for the elementary field (p, with the replacement 
X — s> 0^. From (|T6|) the leading correction for small p is evidently 



SkV^^ = —{ml + Xnx) f ,fj,, , \ ^ k{p) (21) 

4 J {27ry J^^^^^^sXrx 
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which is clearly of order fi^ at the extremum (where ?7i|j + A/jx ~ 0{h))^ showing that 
the extremal energy is independent of the choice of smearing function to this order. 
The actual value of the composite field in the ground state is shifted by the smearing. 
The amount of the shift (extremizing V to leading order in h) is easily seen to be 

where = m\ + SA^jXex is the physical scalar mass (xex is the extremum value of 
the composite field). To understand the origin of this shift we can write 

< 0| y" dr<j){x + r/2)K{r)(j){x - f/2)|0 >=< >^ + / dfK{r)AF{r) (23) 

The second term is just the connected contribution to the expectation value of the 
composite bilocal operator, i.e the propagator for the physical mode smeared with K. 
Going to momentum space and integrating out the energy component, we find 



Thus we see that the shift in < x > at the minimum x = induced by the 

smearing is exactly equal to (0). In general, the composite field can be interpreted 
as equivalent to the square of the usual elementary order parameter as long as the 
smearing is carried out over scales much larger than the inverse mass gap of the 
physical modes. The above argument also makes it clear that the domain of defini- 
tion of the effective potential, for a fixed finite smearing scale p, actually begins at 
X = Xmin > where Xmin is, from the Lehmann representation, at least as large as 
^ -f Tj^-^iP) — / 2 i 1 is (finite) residue of the renormalized propagator at 

V "^i +™-ph 

the physical pole = rUp^. For p small, this end-point can be made as close as we 
wish to zero. 

The upshot of the preceding discussion is simply this: the smearing dependence 
of the effective potential defined with a smeared composite field will be exponen- 
tially small as long as the scale of smearing is kept large compared to the Compton 
wavelength of the physical scalar. Once the field points are sufficiently far apart, the 
expectation of the composite field reverts by clustering to the square of the vev of the 
elementary field, with exponentially small corrections arising from the connected part 
of the bilocal expectation value. This means that ambiguities arising from various 
choices of the bilocal smearing are (a) physically well understood, and (b) reducible 
to an arbitrarily small value. 
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3 Composite Effective Potential for the Higgs- Abelian 
model- Coulomb Gauge 



In this section we shall repeat the calculation of the 1-loop composite effective po- 
tential in a physically more interesting case, that of a spontaneously broken gauge 
symmetry. As we wish to work with gauge-invariant quantities while retaining a 
strict energy interpretation for the potential, we shall use a gauge-invariant bilocal 
operator as our probe of symmetry breaking, but work in a physical gauge where a 
positive Hamiltonian and positive metric Hilbert space obtain. In such a situation 



the standard argument |jT6|] shows that the value of the effective potential at any field 
value is the minimum energy compatible with that value in a physical state (with 
the usual proviso that the effective potential is the convex hull of the perturbatively 
computed Legendre transform). Either Coulomb or axial gauge would be suitable for 
this purpose. The calculation is somewhat simpler in Coulomb gauge, where we can 
maintain the smooth smearing procedure used in Section 2 for the discrete symmetry 
case, so we shall use this gauge henceforth. 

The Lagrangian is a Higgs-abelian model, with an additional fermion (the top 
quark, in the electroweak context) coupled chirally as in the sigma-model: 

£ = |(9^+^eA^)0|^-^F^2^-mV*0-P(0*0)+^^(^-^|754)^-^?j/^(0i+^7502)V' (25) 

where the scalar field = -^{(pi + i(p2) is now complex. The bilocal operator used to 
define a potential will now be taken as 

X(t, x) = J df(P%t,x + r/2)K{f^e'I '^'y^'^^'^'^^-^^^'y^c^it, x - f/2) (26) 

where the c-number current J7'(x, y, r) satisfies 

V^-J{x,y,r) = 5{x-y-f/2)-5{x-y + f/2) (27) 
J{x,y,r) = Vga (28) 

Here <j{x) is a c-number scalar field which is formally identical to the electrostatic 
potential in a dipole field (where J corresponds to an electric dipole field). The choice 
of J ensures that x is a gauge-invariant field. Moreover, the gauge-invariantizing 
phase factor vanishes in Coulomb gauge, by a spatial integration by parts in the 
exponent. Also, the composite field x contains by explicit construction only fields on 
a single time-slice, so the Schrodinger picture essential for the energy interpretation of 



the potential is maintained. Rigorous arguments show |17] that the bilocal operator 
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with a smeared gauge-invariantizing string of this type is a vahd order parameter 
for SSB in this theory, as the expectation value of this bilocal operator remains 
nonzero even when the field points are taken infinitely far apart. Finally, insertion of 
a complete set of states (in Coulomb gauge we have a positive metric Hilbert space) 
leads to positivity of < 0|x|0 > as for the simple theory of Section 2 provided the 
Fourier transform of the smearing function K{r) is positive. In a slight change of 
notation from the preceding section, henceforth bare masses and couplings will be 
unsubscripted, and we no longer display explicitly powers of h. 

As we are dealing with an abelian theory here, the functional integral over the 
gauge degrees of freedom is gaussian, and it will be convenient to perform this inte- 
gration (in Coulomb gauge) immediately. The integration over the fermion field may 
likewise be performed at the outset. The fermion determinant gives a contribution 
to the effective action which is explicitly of order h, and which is a functional of the 
scalar and gauge fields. The gauge field dependence of this functional only contributes 
at order and may therefore be neglected when the gauge field integrations are per- 
formed. (The scalar field dependence of the fermionic determinant will be computed 
and included below.) The remaining functional integral over the scalar field will then 
be evaluated using a saddle-point expansion as in Section 2. In doing the gauge inte- 
grations, it is best to separate the transverse and Coulomb modes, which give rise to 
physically distinct contributions. Thus, defining 

= ie(0*9'^0 - 09'^(/)*) (29) 
D = + 26^*0 (30) 

we find for the integral over the transverse modes 

- i In (Ztr) = -i In y VAS{W ■ l)e^/(-H^^+'^"^^)<i'^ (31) 

= ^1 d''xJiD-'/\5ij - ^)D-'/^Jj + itT\nD (32) 

while the integral over the Coulomb mode yields 

-iln(Ze) = _iin|p^e^/(X-^+2eV<^Mo-JoAo)d^x (33) 

Note that terms quadratic in J are effectively tree-level in powers of Planck's constant, 
and must therefore be included in the effective scalar action when we perform the 
saddle-point evaluation of the effective potential. 
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After integrating out the gauge degrees of freedom, the effective tree-level scalar 
Lagrangian, augmented by a source term for the bilocal x field, becomes 

where = + jK incorporates the source bilinear in the field together with the 
bare mass term into a single nonlocal kernel. Now we expand (p = (po + (p, where 
-Q^{4> = 4>o) = and take the source j{x) and consequently the response (poi^) = 4>o 
constant. One finds 

CtM) ^ 4M + ^(0* ^)M, I ) (36) 

Note that only the Coulomb part of ( P^D contributes in the translationally invari- 
ant limit to the quadratic form A^^, as in this limit Ji —>■ ie{(pQdi(p — (podiCp*) and 
(5,, - = 0. Explicitly, 

= 

g0*2|(/.o V^O-A+2e2|<AoP \^^nL) Q^Q^*\<t>o^ ^ \'V0\ _A+2e2|(/.o|2 

In I r^2\ a^p , I 2 1 J, 12 dl d^p, 2 j,*2 9^ 

-t m ) o^9^*\4>o \(P0\ -A+2e'^\M'^ 902 I (/.Q <^ <P0 -A+2e2|</,o|2 

Specializing to the case P{(p) = \{(p*(pY, and noting that in the translationally 
invariant limit X = I d'^x df(p*{x — r/2)K{r)(p{x + r/2) = |0oP; this simplifies to 
M^ = 

-(2A + e'^Sl?^)<Pl -(□ +rh^ + QXx) + (2A + e'^^)x 



-A+2e2x''^0 ^ -r w/v^; -r v^'v -r <^ -A+2e2x' 

whence the integral over the fluctuation fields (p yields directly 



q2 

lndet(X0) = In det{(n + + 6Ax)^ -(□ + + 6Ax)(4Ax + 26^X3^^^^^)} 

= Indetfn + + 6Ax) + IndetiD + rh^ + 2Xx - 2e^X—. — ^-^} 

—A + 2e^x 

= lndet(n + + 6Ax) + lndet(n + 2e\) - lndet(-A + 2e\) 
—A + 2e^Y 1 

+ lndet{l + ( — -) ^(m2 + 2Ax)} + constant (37) 

—A □ + 2e^x 
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Adding in the one-loop contributions from the gauge integrations we find 



|trln(n + 2e2x)+2^1 



^gai°+scaiar = ^-tr lu (□ + 2e2x) + 2-tr lu (□ + + QA^) 



+ + ( _^ -- ) ^^^^^^ (m^ + 2Ax)} (38) 

The unpleasant term lndet(— A + 2e^x)) which contains a UV divergence of the form 
/ dpo not present in the counterterms, has cancelled, and the remaining terms are 
easily seen to correspond in the broken symmetry phase to three massive gauge vector 
modes and a single massive scalar mode. The complicated and peculiar last term is 
a remnant of the long-range Coulomb interaction, and contains divergences which 
will be removed by the counterterms of the theory, as we shall see shortly. In these 
one-loop contributions, the source augmented mass rh contains the source jo of the 
composite field determined at tree level from -g^\<j,o = for 0o constant, so 

Jo = (-m2 + F'(x)) = -K + 2Ax) (39) 
m2 ^ ^2 ^ ^^j^ = m'^(^i_x)~ 2XxK (40) 

Integrating out the fermion field gives a similar contribution (except for the charac- 
teristic change of sign) 

- V,l;^Z = trln (t^ - gyi<l>, + ^75^2)) (41) 
The full effective potential through 1-loop is thus given by 

V{x) = m\ + + Vi^;/3caia. + V^fL-^oT (42) 

This expression is as yet unrenormalized- we must rewrite the bare parameters m. A, Qy 
in terms of renormalized parameters and rescale the fields appropriately, at which 
point the divergences in the one-loop contributions will be seen to cancel completely. 
This calculation will be performed in the next section. 



4 Renormalized Composite Effective Potential for 
the Higgs-Abelian Model 

To renormalize the result ( ^^ obtained above, one may integrate out the energy 
component po in momentum space and then dimensionally regularize the resulting 
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purely spatial integrals (which are then carried out in 3-e dimensions). For example 
the contribution from the massive gauge vector loop was found to be 



Va 



-ii-trln(n + 2e^x) 



—i- 
2 

const + 
3-e 



(2^ 
3 



In (— + 2e^x) Minkowski 



ln(l+ 2-'^ 



d^- 



{2-K 



pi + p^ 



Euclidean 



(27r) 



3-e ' 



]P + 2e'^x - \p\] 



(43) 



where the final spatial integral, together with the prefactor {g^^ — 1) counting spatial 
gauge modes, has been dimensionally continued. The divergent part of Va can be 
separated off as follows: 



Va 

VaMv 



3-e 
2 

3-e 



d^ ^p 
(27r)3-^^ 



4 2 



p- + 2e^X-|p1-^ + 2(|^-^,_^^,^3/.. 



+ Va, 



div 



4 
3 



16n' 



(27r)3-^(|p]2 + ^2)3/2 



(44) 



and 7] is an infrared cutoff which will 
^) has all the subtractions needed to 



Here the divergent piece is defined in (|T/ 

be set to zero at the end. The first term in ^ ^ 

ensure that it remains finite as e — so we may define a MS-renormalized gauge 
boson contribution (obtained as usual by subtracting the part proportional to in 
VA,div) as 



3 I- 



d^p 
(2^ 



jp + 2e^x -\V\- + 



\p\ 2(p2 + r/2)3/2 



lOTT^ 3 

(45) 

Similar manipulations can be used to extract the singular parts of the massive 
Higgs, residual Coulomb, and fermion contributions to the 1-loop composite effective 
potential. One obtains for the massive scalar loop contribution 



V, 



B 



P 



(277)3-^ 



f + AXx + (m2 + 2Ax)(l - K{p)) - \p\) 



647r2 



(46) 
(47) 
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with 

Vb ren 



+ 



(27r)3 
1 (m2 + 6\xf 



f + 4:Xx+ (m2 + 2Ax)(l - k{p)) - \p\ 



m 



^ (m^ + 6Ax)^ln'^ 



8 ^2)3/2 ^ 64^2 . 

The residual Coulomb contribution also yields a divergent contribution: 



(m' + 6Ax) 
(48) 



Vc 



1 r "p 

2 J (27r)3-^ 



e\{m^ + 2Ax) A 



1 



p2 



647r2 



(m2 + 2Ax)'A, 



(49) 
(50) 



This term couples the Higgs and gauge sectors and the renormalized contribution is 
rather complicated. One finds 



l + („2 + 2Ax)^-^^^-l) 



p2 



\ 



p2 



+ 



1 



1 + (m2 + 2Ax) ^ 5^^^ - 1 - 



p2 



2p2 

+ 2Ax 
2(p2 + 772) 



8(p^ + 7?2^2 



)} 



e2x(m2 + 2Ax)(ln(\)-2 + 21n2) 



+ 



167r2 
1 



m2 + 2Ax)'(ln(\)-l + 21n2) 



647r2 A* 
Finally, the 1 loop fermion contribution gives 

d^'^p 



(51) 



D 



-2 



(27r)3-^ 

= VD,ren + 7^«?y A, 



47r2^ 



(52) 
(53) 



where 



^Aren = "2 / ^ ( ^^T^ - |p1 - 



4,, 2 



2\p\ 8 (p2 + 772)3/2. 



47r2 



(54) 
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The 1-loop renormalization of the theory (mass, couphng and field rescaling) im- 
plies the following cutoff dependence for the coefficients of x and in the tree contri- 
bution to the potential (individual counterterm contributions are listed in Appendix 
A): 



m^(l + _l^(-e|, + 4A«)) (55) 



A, 
167r2 

A - \R + ^^{-2el\R+lQ\l + ?,e\-Agl) (56) 



When these replacements are made in the unrenormalized potential (|4^), one finds 
that all divergences cancel and we are left with an explicitly finite result for the 
composite potential (now in terms of MS-renormalized fields, masses and couplings): 

VrcM = + ^RX^ + VA,ren(x) + VB,rcn(x) + Vc,ren(x) + VD,rcn(x) (57) 

with the renormalized potentials given explicitly in ( ^ , |i8| , |51| , ^4] ) with the mass and 
couplings therein interpreted as renormalized ones. 



5 Analytic Formulae for Sharp Momentum- space 
Smearing 

As mentioned previously, any choice of smearing function K (r) which adequately 
suppresses short- distance contributions to {x (x)) could be used in defining the com- 
posite operator. To the extent that different K (r) isolate only the longest- distance 
modes relevant for SSB they should give similar values for {x{x)) and thus express 
the same physics. In particular, they should unambiguously signal the location and 
depth of symmetry-breaking extrema of the effective potential. One simple choice 
which permits closed-form expressions for the dimensionally-regulated momentum 
integrals is the step function, K (p) = 9 {p — \p\) . With this choice, the i^T-dependent 
integrals in Vb and Vc can be written generically (to (9 (e°)) as 

= I d'-'pf{0)+ I d'p (/(l)-/(0)) (58) 
J J\p\<p 

The first integral contains the UV divergence and can be evaluated in the usual 
manner. It is independent of the smearing radius p. The second integral is manifestly 
UV finite, contains all of the dependence on p, and is also straightforward to evaluate. 
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Thus we see that at one loop the effective potential for the bifocal composite operator 
naturally decomposes into a piece which is the elementary field effective potential 
(written in terms of % = 0*0) and a piece dependent on the smearing scale which 
introduces no new UV divergences, unlik;e the local composite operator 0^(a;). 
The explicit expressions for terms of the effective potential are: 



^a(x) = 
^b(x) = 

Vb,2{x) = 



SlTly 

64^ 



-A, + In ^ - - 
/i^ 6 



Vb,i (x) + Vb,2 ix) 



647r2 
1 

64^ ' 



647r2 ' 



1 rnt 3 



(59) 

(60) 
(61) 



2p^p2 + {2p' + 4Ax) - (4Ax)' In ^ — 

^P^f^^H (2P^ + ml) - ml In ^'^^P"^"^"^ 



Vcix) = K^,i(x) + K;,2(x) 



Vb,i(x) = g^|-("^s-"^v)^(A6 + ^j -2mgmy(mg-my)=^ + 



^0,2 (x) = ^ ^P\/^l + H + 2p') - In 



In^r- 



(62) 
(63) 

(64) 



)- 



647r^ j^' V « ' V y ' y « j 

^^j-y/m^ + p^^Jm"^ + p2 ^77^2 j^jy^^j^ 2p^) - mgmv 

( 2 2\2, v/p^ + ^^ + Vp^ + < | 

I m — my ] In > 

^ " ' mg + mv J 



[m] + ml) - 



1 "^f 3 



(65) 
(66) 



where = m^ + 2Ax, = 2e^X) m\ — m^ + 6Ax and = 2g'^x. 

As noted previously, for p — > the p-dependent terms go smoothly to zero. For 
p — > oo new quadratic divergences arise as the point-split composite operator ap- 
proaches the local composite operator. 
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6 Renormalization Group Improvement and Vac- 
uum Stability Bounds 



The composite operator effective potential formalism developed above provides a 
gauge-invariant framework for obtaining the information about symmetry breaking 
which is partially obscured by gauge dependence and UV problems in other treat- 
ments. While the value of the elementary-field effective potential is gauge invariant at 
any of its local extrema, the field value at which the extremum occurs (i.e. the expec- 
tation value of the elementary scalar field in the corresponding phase of the theory) is 
not. Thus, associating the gauge dependent expectation value of the elementary field 
at which some feature of the effective potential occurs with gauge invariant physical 
quantities is a practice of questionable validity. The expectation value of the compos- 
ite operator is by construction free of these gauge ambiguities. Because the effective 
potential of the composite operator has an energy interpretation, the value of the 
effective potential at any given value of the composite order parameter corresponds 
to the minimum physical energy density compatible with that value. Consequently, 
the global minimum of this effective potential should correspond to the true vacuum 
of the theory and the value of x at which it occurs should also have physical meaning. 

One context in which such a treatment is useful is in the formulation of a lower 
bound on the Standard Model Higgs mass from vacuum stability considerations. If 
the electroweak vacuum is assumed to be stable (rather than merely metastable) 
then it must be the global minimum of the effective potential (in both the case 
of the composite operator effective potential and of the elementary field effective 
potential). In practice, however, we expect that the Standard Model is the low- 
energy effective theory of a more fundamental high-scale theory and so will only 
be an accurate description of physics up to some energy scale (the 'new physics' 
scale, perhaps characterized by the masses of new heavy particles). Thus, it is only 
consistent to demand that the electroweak vacuum is the global minimum up to the 
scale at which the effective theory breaks down, beyond which the model is simply 
no longer accurate. 

For the elementary field effective potential, the statement that the value of the 
effective potential is larger than the value at the electroweak minimum up to some 
scale A is a gauge- dependent statement @]. As such it is an unsatisfactory criterion for 
defining a 'new physics' scale. For the composite operator effective potential defined 
and evaluated in this paper there is no gauge dependence and we may interpret the 
composite operator expectation value as a physically meaningful energy scale. 
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Since the Standard Model may be valid up to very high energies it will be necessary 
to study the effective potential for field expectation values much larger than the elec- 
troweak scale. However, it is well-known that in order to study the effective potential 
at large field values, the usual perturbative loop expansion for the effective potential 
may be inadequate. For the elementary field effective potential, the loop expansion 
generates terms generically of form In (where g is any of the couplings) . If 
the couplings g and the renormalization scale fi are given values characteristic of the 
electroweak scale (e.g. /i = Mz and the couplings given their values at the Z-scale), 
considering ^ /i will invalidate the perturbative expansion. The range of validity 
of the approximation may be improved by utilizing the renormalization group (RG) 



improved effective potential [1^, |19[. Similar considerations apply to the composite 
operator effective potential. 

The full (all-orders) effective potential is independent of /i. This independence 
can be expressed as a first-order differential equation 

Ai^Kff (x,m,e, A,5(j^,/i) = (67) 

This equation can be solved and, combined with dimensional analysis, yields the 
RG-improved effective potential @, |l^ 



Vcs (s^Xi, 9i, /i) = [C (s)]^ Kff (xi, g (s, gi) , m {s, rrii) , /x) (68) 

where 

1 Li dZy , , 

Note that 7,^ in (69) is really one half of the anomalous dimension of the gauge- 
invariant composite operator x- In Coulomb gauge (and only in Coulomb gauge) the 
gauge-invariantizing string factor in x is equal to unity and the operator therefore 
renormalizes as the square of the elementary field in this gauge. Accordingly 7^ is 
just the usual anomalous dimension of the elementary scalar in Coulomb gauge. In 
other gauges, there will be a divergent contribution from a nontrivial string factor 
which compensates for the gauge- variance of 7^. The scale factor in (68) is just 

r-log s \ 



C (s) = exp 



7<A i.^) + 1 



dx 



(70) 



S2 

g are solutions to the equations 



Here = ^ and g represents the set of couplings gy, A, and e. The running couplings 



'^^ = M9{s)) = , A(0) = A. (71) 
as 1 + 70 
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and analogous for the other couphngs. The gi are the values for the couplings at 
the initial scale Xi? which we will take to be around the electroweak scale. The (3 
functions can be computed (in a loop expansion) from the MS counterterms of the 
theory and are listed in Appendix A. Whereas the unimproved effective potential was 

2 

trustworthy only for x such that ^^gfs'^ <^ 1, the RG improved effective potential 
will be reliable as long as the running couplings remain small. 

It has been demonstrated (in the context of the elementary field effective potential) 
that the n-loop effective potential improved using n + 1 loop (3 and 7 functions 
resums the n^^^-to-leading logs (^D], |2T|]. Since we are demonstrating the utility of a 
calculational tool rather than pursuing a precise numerical result, we will be satisfied 
to sum the leading logs only. It is thus sufficient to consider the one-loop effective 
potential with the tree-level piece run with one-loop (3 functions. In the large-field 
{x ^ ^^) limit the RG improved effective potential can be written 

Ves (s^Xi, 9i, /Lt) = Aeff (Xi, 9 (s, Qi) , /i) [xiC (■«)] (72) 

Acfr(s) = A(s) + AA(^i,Xi,/i) (73) 

We will demonstrate the calculation of the vacuum stability bound using the RG- 
improved composite operator effective potential in the context of a toy model, the 
abelian Higgs model coupled to a fermion. The issue of vacuum stability arises in 
qualitatively the same way in the abelian Higgs+fermion model as in the Standard 
Model, as a result of the large negative contribution of top quark loops to Px. For 
some electroweak-scale boundary condition on A(s), X{sew) = the top quark 
term in f3x drives Aefr(s) negative for sufficiently large s. For large s the tree-level 
RG-improved effective potential Ves{xi^^) ~ Aeff(s)XiC('5)'^ then falls rapidly below 
VesisEw) and the electroweak vacuum is unstable. We define Syi as the value of s 
at which at which the RG improved effective potential equals zero. Since ({s) > 0, 
Vesisvi) = is equivalent to Xesisyi) = 0. For the electroweak vacuum to remain 
the global minimum of the theory, the 'new physics' must enter before syi (however, 
5S| for qualification). For a fixed Syi there corresponds a minimum Aj below 



see 



which the electroweak vacuum is destabilized too early. This then translates to a 
lower bound on the Higgs pole mass. 

7 Numerical Results 

We present here some numerical results to demonstrate the qualitative features of the 
composite operator effective potential and its use in the study of vacuum stability 
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Figure 1: Effective potentials in the gauge for several values of the gauge parameter. 
The arrows indicate the locations of the minima. 



bounds. A set of initial values of parameters has been chosen (ef = 0.15, Qy = 
0.5, V = 246 GeV, /i = f ) to resemble the Standard Model, but these plots should not 
be construed as a serious attempt to calculate Standard Model quantities. Figures 
|T] and illustrate the gauge dependence of the elementary field effective potential 
(prior to RG improvement) for Aj = 0.2. Fig. |I] shows the effective potential in the 
gauge for several values of the gauge parameter and explicitly indicates the shift 
in the location of the 'electroweak' minimum as a function of the gauge parameter 
^e^. In Fig. 0the one-loop corrections are isolated to highlight the gauge dependence. 
In Fig. ^ we plot the one-loop composite operator effective potential with p- space 
smearing function K{p) = 6{p— \p\) for several values of p. Since the mass parameter 
m, which introduces a typical energy scale for the theory, has been taken to be 
110 GeV, we expect that for p much smaller than this the composite operator will 
be insensitive to all but SSB effects, which are purely infrared. For p much larger 
than this the composite operator will detect the shorter-wavelength fluctuations not 
associated with SSB. As is shown in the figure, curves with p of order the mass 
parameter are very close to the p = curve. Only for p ~ 1000 GeV does the 
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Figure 2: One-loop corrections to the effective potentials in the gauge for several values 
of the gauge parameter. 
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Figure 3: Effective potentials for the composite operator using a sharp p-space cutoff smear- 
ing function. For cutoff momenta of order the mass parameter (m = 110 GeV) or smaller, 
the curves are essentially indistinguishable. 

shape of the composite operator effective potential deviate significantly from the 
p = curve, so we are confident in the ability of the composite operator to provide 
a reliable separation of the physical scales of the problem and to isolate only the 
SSB effects. Fig. § shows the results of a vacuum stability bound calculation using 
the composite operator effective potential with p = 0. For a choice of Syi apply the 
boundary condition Xes{svi) = X{svi) + AX = and use the RG equation to run A(s) 
down to s = 0(x = Xi) ("W'e neglect the mass parameter in our expression for AA). 
The resulting Xi is the minimum value for which the RG-improved effective potential 
remains positive up to scale syi- In a complete treatment this lower bound on Xi 
could be converted to a lower bound on the Higgs pole mass. We have performed 
a similar calculation for the elementary field effective potential in Landau gauge 
and find the results virtually indistinguishable. For the case studied here, with weak 
gauge couplings, this is simply a reflection of the small quantitative contribution from 
the gauge sector, and not a statement of identity of Landau and Coulomb effective 
potentials. This suggests that in the Standard Model as well the Landau gauge 
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Figure 4: Aj^min vs. syi for the composite operator effective potential with p = 0. The 
corresponding curve obtained using the elementary-field effective potential in Landau gauge 
is virtually indistinguishable. In this plot, ef = 0.15, gy = 0.5, v = 246GeV,;U^ = Xi = ■ 

results may be similar to those obtained using a gauge-invariant method. However, 
due to numerical differences in the j3 functions, the effect of QCD in the running 
of Qy, and of course the different gauge group of the Standard Model, no immediate 
conclusions can be drawn in that context. In Fig. |^ we plot the one- loop correction to 
the elementary-field effective potential in Landau gauge and the composite operator 
effective potential with (p = 0). The fact that the curves are so close in the region 
in which we choose initial scale Xi explains why the corresponding AA and thus the 
Ai,min curves are so similar for the two effective potentials. 

8 Composite Operator Effective Potentials for Non- 
abelian Gauge Theories 

In an abelian gauge theory, the use of Coulomb gauge was facilitated by the ob- 
servation that a smeared string corresponding to a dipole field (cf. (25-27)) yields a 
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Figure 5: One- loop corrections to the elementary field effective potential in Landau gauge 
and composite operator effective potential (for cutoff momentum zero), for = 0.15, 
= 0.5, A = 0.2, ji = V. Both curves are without RG improvement. 
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gauge-invariantizing factor for a bilocal operator which automatically reduces to unity 
in Coulomb gauge. The need for path ordering of such factors in the nonabelian case 
means that this procedure fails and we are forced to seek a more convenient physical 
gauge. At first sight, axial gauge would seem to fit the bill- in this gauge, there is a 
positive definite state space and a well defined Hamiltonian, and the ordinary straight 
line string factor can be used to define a gauge-invariant composite field in the usual 
way 

X,D{t,x)= f dl^* (t,x+^-z] K{l)Ve'I'^^'^^^'>'(P* (t,x-^-z] (74) 



Here the scalar fields have been fixed to lie along an arbitrarily chosen line z. The 
path ordering along the one- dimensional (spacelike) path is unambiguous. If the 
calculation is then performed in the A-z = spacelike axial gauge the string operator 
again reduces to unity. 

The use of axial gauge in this fashion leads however to a number of interesting 
delicate points, which divide into two categories- perturbative ultraviolet problems 
in defining the effective potential, and nonperturbative infrared problems associated 
with the loss of long range order when such string operators are employed as 
order parameters in a broken gauge theory. 

First, we discuss the problems arising with the definition of a finite composite po- 
tential in perturbation theory. These problems are not specific to gauge theory - aris- 
ing simply from the UV singularities of one-dimensionally smeared bilocal operators- 
and can be illustrated in the simple scalar model of Section 2. We found there that 
the one loop composite effective potential yields a term 

/ ^(Vp^ + 2Xnx + {ml + XnX){l - K{p}) - \p\) (75) 

where K{p) is the Fourier transform of the coordinate space bilocal smearing function 
K{r). It is crucial that the terms proportional to powers of K in the expansion of this 
expression in powers of x iT'Ot introduce additional ultraviolet divergences, as there are 
no counterterms available in the Lagrangian to absorb them. If we simply point split 
the fields (along the z direction, say) by taking K{r) = ^6{x)S{y){d{x — 1/p) + 6{x + 
1/p)), corresponding to = cos y, we find that the one- loop potential expanded in 
X is UV finite at 0{x) but logarithmically divergent at (9(x^)- Namely: 

d p I < oo (76) 

]? + m% + 3Xrx 



I 



d Pt^ 9 r x./o ~ - / d^PT^ ^ ^ + finite (77) 
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~ log divergent 



(78) 



The origin of the divergence is not hard to find- at order x^, the momentum space 
double insertion of the point split operator (f){r+^z)(f){r — ^z) leads to a UV divergence 
when the coordinate space integration defining the Fourier transform brings both pairs 
of field operators simultaneously together. Clearly, holding the two fields at a fixed 
separation is a prescription for trouble. On the other hand, taking a smooth smearing 
function such as = e"^^''^''^ (in analogy to the 3 dimensional Gaussian smearing of 
Sections 2,3), leads to an uncompensated divergence at 0(x), with order and higher 

terms finite. The linear divergence in x / d^p—^=^^i^ is also readily understood 

in operator terms. Here the one-dimensional smearing integral leads directly to a 
linear divergence due to the quadratic short- distance divergence of a product of two 
scalar field operators approaching the same point. 

The solution to these ultraviolet problems is actually quite simple- one need only 
choose a smearing function which vanishes at the origin in coordinate space. This 
removes the divergence from the one dimensional integral in the region where the 
two field points approach one another in a single insertion of the composite field 
while removing the logarithmic divergence at order by smearing out the split field 
locations. For example, a suitably normalized smearing function is 

K{r) = -^p^6{x)6{y)z^e"p''' (79) 

As for the three dimensional smearings of Sections 2-4, one should take the smearing 
parameter p small, corresponding to widely separated peaks in K{r). In the opposite 
limit of p large, the UV divergences will reappear, distorting the shape of the effective 
potential. Of course, the requirement that K{0) = now means that the Fourier 
transform K is not positive for all p, so there is no rigorous argument that the natural 
domain of x is restricted to the positive real axis. Still, for small p in the symmetry- 
broken phase, we expect that the minimum of the effective potential will be found at 
a safely positive location, so this property is perhaps not terribly important. 

A further difficulty which one encounters in employing xid as an order parameter 
is specific to gauge theories with nontrivial topological structure. The decorrelating 
effects induced by instantons on bilocal operators with a gauge-invariantizing string 
was first pointed out by Frohlich et. al. pGf. The existence of Gribov copies can 



also lead to a destruction of long range order, again as a result of large nonperturbative 
field configurations. All such effects are presumably of order e~^^^ and therefore 
not visible in a standard perturbative loop expansion. From a practical point of 
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view, in the weakly coupled electroweak case instanton effects (unenhanced by large 
combinatoric prefactors as is potentially the case in multiparticle production at high 
energy) are extremely small, so these effects are clearly ignorable. In strongly coupled 
theories- when studying dynamical symmetry breaking, for example- the problem 
returns and caution will be required when taking the limit p ^ in an axial type 
gauge. We should remind the reader that for the Coulomb gauge smeared bilocal 
operator, Kennedy and King have shown that this limit leads to a nonvanishing 
order parameter in the symmetry broken phase. 



9 Conclusions 

The computation of lower bounds on the Higgs mass from vacuum stability constraints 
using the elementary field effective potential results results in unphysical gauge de- 
pendence in such bounds. We have formulated an effective potential in terms of gauge 
invariant composite operators which avoids the ultraviolet problems of local compos- 
ite operators while retaining an energy interpretation. We have shown how this could 
be used to calculate the vacuum stability bound on the Higgs mass in the context of 
a toy model, the abelian Higgs model coupled to a fermion, and for the set of param- 
eters chosen (corresponding to a weakly coupled gauge sector) find that the results 
are quantitatively close to those obtained from the elementary field effective potential 
in Landau gauge. While the extension to a nonabelian gauge theory using a related 
composite operator in axial gauge introduces additional subtleties, we see no serious 
obstacles to the complete calculation. Thus, the tool might be extended to models of 
phenomenological interest, such as the Standard Model or its extensions. Composite 
effective potentials are essential in the study of dynamical symmetry breaking, so the 
issues discussed in this paper should prove of value in that context as well. 
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Appendix A: MS renormalization 

MS counterterms 



5Z^ 
5Z^ 



5Z, 



9y 



5Ze 

5Za 
5Zt 



(47r)^J 



-10A + 66^-3^-4^2 ^4^ 



A 

-A. 



(-4A + 36^ - 



— -e 



L(47r)^ 
A, 



L(47r) 
-A/ 



- L(47r)^ 



(47r)^J 



-Ae 

(47r)^ 



/3 functions 



The relevant one-loop MS (3 and 7 functions for the theory are: 
/3a = ^(20A2 + 6e^-8^^-12e2A + 8Ayj) 

= 16^ (""' + ^^) 
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